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NOMENCLATURE 


specific  heat,  BTU/lbm-°F. 
least-squares  error,  Equation  (5),  p.  21 

effective  thermal  conductivity  of  composite,  BTU/(hr-°F-f t) 
temperature,  °F. 
time,  hour  or  second 
volume  fraction 

orthogonal  coordinate;  x  refers  to  the  heat  flow  direction. 
Figure  3 

vmbols 

effective  thermal  diffusivity,  ft2/sec 
density,  lbm/ft3 


Pts 

refers  to  fiber 
refers  to  resin  matrix 
along  fiber  axis 


transverse  to  fiber  axis 
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INTRODUCTION 

Among  che  composite  materials,  those  with  fibers  in  a  matrix  are 
the  most  popular  in  aerospace  structures.  There  are  several  reasons  for 
this  popularity:  one  is  the  high  strength  and  modulus  of  fiber  which 
contributes  significantly  to  the  overall  strength  of  the  resul  g  com¬ 
posite  materials.  Next  is  the  low  thermal  expansion  of  these  >rs 
which  manifests  itself  in  the  resulting  dimensional  stability  ruc- 

tures  so  formed.  Furthermore,  with  the  advent  of  pre-impregnated 
(prepreg)  tapes,  formability  into  various  shapes  is  greatly  enhanced. 
Thus,  the  increased  use  of  fibrous  composites  is  coupled  with  lower 
cost;  their  expanding  applications  from  aerospace  to  other  adaptations 
are  an  Indication  of  their  potential  role  as  a  candidate  material  for 
many  applications.  Composites,  as  a  new  class  of  materials,  have 
emerged  to  the  full  status  of  an  engineering  material. 

Research  efforts  in  the  past  have  chiefly  focused  in  the  mechanical 
aspects  of  composites— manufacturing  processes  and  attendant  strength, 
etc.  Recent  potential  applications  in  which  thermal  considerations 
have  become  a  crucial  factor  point  the  need  for  better  tools  In 
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analyzing  Che  thermal  performance  of  composite  materials  in  general  and 
fibrous  composites  in  particular  under  various  heat  loading  conditions. 

As  in  the  case  of  research  into  the  mechanical  aspects  of  composite 
materials,  Investigations  with  regard  to  the  thermal  aspects  proceeded 
along  several  avenues.  In  specific  industrial  applications,  once  a  com¬ 
posite  structure  had  been  designed  or  fabricated,  basic  thermal  perform¬ 
ance  parameters — such  as  Che  apparent  thermal  conductivity  and  a  few 
others — were  usually  obtained  in  situ.  The  resulting  data  were  only 
specifically  addressed  to  that  composite  material  studied.  Data  ob¬ 
tained  on  such  an  ad  hoc  basis,  scattered  in  various  industrial  labora¬ 
tories,  were  not  amenable  to  overall  coordination  in  terms  of  basic  cor¬ 
relations  between  the  geometrical  dispersion  pattern  of  fibers  in  matrix, 
the  constituent  conductivities,  and  the  apparent  thermal  conductivity. 

Parallel  to  the  above-mentioned  type  of  measurements  for  their 
specific  needs,  basic  investigations  of  heat  conduction  in  composite 
media  have  been  proceeding  in  the  past  decade.  Wi  hln  this  framework 
of  basic  study,  two  approaches  were  pursued:  (1)  one  was  to  devise 
methods  of  calculation  of  the  temperature  distribution  in  composites 
with  simplified  dispersion  patterns  and  for  shapes  of  simplified  geome¬ 
tries.  Inherent  in  this  approach  of  analysis  was  that  the  thermal 
physical  properties  of  the  constituents  in  a  composite  medium  were  known 
a  priori.  This  type  of  investigation  was  exemplified  by  the  work  of 
Horvay,  et  al.  [l],  who  addressed  the  problem  of  transient  heat  flow  in 
the  plane  of  laminates  for  a  layered  composite.  Later  on,  Manaker  and 
Horvay  [2]  extended  the  original  analysis  to  problems  with  various 
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heating  boundary  conditions.  These  were  but  a  small  sample  of  a  good 
number  of  publications  along  this  approach.  As  important  and  elegant 
were  these  investigations,  the  needs  of  design  engineering  in  estimating 
thermal  conduction  in  composites  were  not  fulfilled;  for  want  of  the 
thermal  physical  properties  of  the  constituents;  for  lack  of  a  correla¬ 
tion  between  the  input  parameters  (the  thermal  properties  of  the  consti¬ 
tuents  and  the  patterns  of  dispersion) ;  and  for  need  of  the  overall  ther¬ 
mal  performance,  (ii)  The  second  approach  was  to  model  a  composite  medi¬ 
um  as  a  single-phase  body  which  has  thermal  characteristics,  in  particu¬ 
lar  thermal  conductivities,  reflective  of  those  of  the  ingredients.  In 
other  words,  the  principal  thermal  properties  were  expressed  by  an  effec¬ 
tive  thermal  conductivity,  which  may  be  directionally  dependent,  and  an 
effective  thermal  inertia.  These  effective  values  were  manifestations  of 
Chose  of  the  constituent  materials,  in  the  case  of  a  fibrous  material, 
the  fibers  and  the  matrix  material.  The  relationships  between  the  con¬ 
stituents  and  an  effective  medium  as  a  model  replacement  were  the  main 
focus  of  many  research  efforts.  The  establishment  of  such  a  relation¬ 
ship  was  a  principal  objective  which  resulted  In  a  publication  by  Han 
and  Cosner  [3]. 

Since  the  previously  described  efforts  of  both  types  were  predi¬ 
cated  on  a  foreknowledge  of  the  basic  thermal  properties  of  the  consti¬ 
tuents,  a  natural  question  arose  as  to  the  availability  of  the  informa¬ 
tion,  in  particular,  relative  to  carbon  fibers  and  resin  as  a  matrix 
material.  Data  of  this  sort  were  meager  in  the  literature.  Consequent¬ 
ly,  a  purpose  of  this  experimental  investigation  within  the  framework  of 
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a  larger  effort  to  develop  heat  transfer  methods  for  composite  material 
was  to  provide  Information  needed  to  bridge  the  gap.  A  secondary  pur¬ 
pose  of  this  experimental  Investigation  was  to  evolve  a  reliable,  inex¬ 
pensive  experimental  procedure — by  way  of  demonstration  with  carbon/- 
epoxy  composites  so  that  other  composites  may  be  similarly  treated. 

With  these  two  objectives  in  mind,  the  geometrical  dispersion 
pattern  of  fibers  in  a  binding  matrix  was  made  as  simple  as  possible  and 
more  importantly  to  conform  with  the  dispersion  pattern  analyzed  in  [3]. 
Hence,  a  carbon/epoxy  composite  with  a  uni-directional  fiber  orientation 
was  used  in  this  study. 
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II 


BASIC  CONSIDERATIONS 


In  an  isotropic  medium,  heat  conduction  is  governed  by  three  phy¬ 
sical  properties  of  the  medium:  thermal  conductivity  k,  density  p,  and 
specific  heat  C.  The  last  two  are  invariably  joined  together  as  a  pro¬ 
duct  term.  From  these  three  parameters,  the  thermal  diffusivity  o  is 
defined  as: 


a  -  k/(pC) 

For  steady-state  conduction,  the  thermal  conductivity  k  alone  determines 
the  entire  temperature  distribution  throughout  the  body;  whereas  for 
transient  heat  flows,  the  thermal  diffusivity  coefficient  a  becomes 
equally  important.  Conceptually,  it  denotes  the  ability  of  diffusing, 
or  passing  on  by  virtue  of  conduction,  that  energy  originating  at  the 
input  site  in  a  material  body  in  question. 

It  is  of  Importance  to  note  that  the  physical  properties  p  and  C 
are  scalars  and  the  thermal  conductivity  k  may  be  directionally  depend¬ 
ent,  in  which  case  the  term  orthotropic  is  used.  For  composites  consis¬ 
ting  of  uni-directional  fibers  in  a  body  of  the  binding  matrix,  transfer 
of  heat  encounters  different  resistances  depending  on  the  direction  of 
heat  propagation.  Hence,  in  composites,  there  are,  generally  speaking. 
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three  effective  thermal  conductivities:  one  for  heat  conduction  parallel 
to  the  fiber  axis,  and  the  other  two  for  transverse  heat  flows  in  the  re¬ 
maining  directions.  Assuming  (in  fact,  true)  that  carbon  fibers  are 
more  conductive  than  the  matrix  materials,  conduction  of  heat  in  the 
direction  of  the  fiber  axis  therefore  constitutes  a  series  of  parallel 
flow  circuits  with  the  flow  along  the  fibers  as  a  short  circuit.  High 
effective  conductivities  are  a  consequence.  For  heat  conduction  trans¬ 
verse  to  the  fiber  axis,  the  high-conductivity  fibers  become  local  in¬ 
clusions  embedded  in  the  low-conductivity  matrix  material;  therefore,  a 
lower  effective  conductivity  is  expected.  Aside  from  these  effects, 
carbon  fibers  and,  for  that  matter,  most  synthetic  fibers,  exhibit  dif¬ 
ferent  conductivities  when  heat  conduction  occurs  parallel  versus  trans¬ 
verse  to  them.  That  is,  due  to  Che  manner  in  which  the  poly¬ 
mers  are  aligned  with  the  fiber  axis,  a  higher  conductivity  results 
when  heat  is  transferred  parallel  to  the  fiber  rather  than  normal  to  it. 

A  cross-sectional  view  of  a  sample  of  the  carbon/epoxy  composite 
used  in  this  study  is  shown  in  Figure  1  indicating  the  discrete  fibers 
(more  properly,  bundles  of  fibers)  in  a  resin  matrix.  The  two  groups  of 
fibers  are  separated  by  a  wide  channel  consisting  of  the  matrix  material 
only;  the  appearance  of  a  canal- like  structure  is  a  result  of  two  prepreg 
tapes  interfacing  each  other.  Within  each  tape,  there  is,  however,  an  or¬ 
derly  identifiable,  geometrical  pattern.  The  cross-sectional  view  of  Fi¬ 
gure  1  was  obtained  by  mounting,  polishing  and  then  microscopically  photo¬ 
graphing  an  end  view  of  the  composite  specimen.  A  further  enlargment  of  the 
microscopic  image,  2100  times  the  original  size,  was  then  used  in  order 


6 


Co  determine  Che  volume  ratio  of  Che  composite,  i.e.,  Che  rado  of  Che 
fiber  volume  Co  Che  cocal  volume;  details  were  documented  in  Appendix  B. 

1.  EXPERIMENT  RATIONALE 

In  principle,  the  thermal  characteristic  of  a  composite  consis¬ 
ting  of  uni-directional  fibers  in  a  matrix,  are  completely  defined  once 
the  thermal  properties  of  the  constituents  and  the  geometrical  disper¬ 
sion  pattern  are  known.  Indeed,  such  was  the  basic  tenet  in  a  previous 
effort  [3].  However,  direct  measurements  of  the  properties  of  the 
elementary  materials  entail  a  great  deal  of  difficulties,  especially 
for  the  fibers  because  of  the  directional  dependence.  In  this  experi¬ 
mental  investigation,  a  transient  heat  conduction  technique  was  used, 
first  to  determine  the  effective  thermal  diffusivities  in  the  three 
principal  directions,  and  then  to  synthesize  the  data  obtained  in  order 
to  delineate  the  thermal-physical  properties  of  the  constituents. 

Conceptually,  the  experiment  consists  of  imparting  a  constant 
heat  flux  on  one  surface  of  a  composite  specimen  block  while  the  surface 
opposite  to  the  heating  surface  was  kept  at  the  initial  tmperature  of 
its  specimen  block.  With  the  remaining  peripheral  surfaces  Insulated 
against  heat  losses,  one-dimensional  heat  conduction  was  thereby 
achieved,  for  which  the  governing  differential  equation  was: 

pC(3T/3t)  -  k(32T/3x2)  (1) 


where : 


8 


T  ■  local  temperature  of  the  composite  body 
t  ■  time 

x  "  coordinate  of  the  heat  flow  direction 
k  »  effective  conductivity  for  heat  flow  In  the  x  direction 
pC  “  effective  thermal  inertia  of  the  composite  material 

The  temperature  T  in  Equation  (1)  represents  an  average  value  of  the 
fibers  and  matrix  at  the  location  concerned. 

Since  the  product  (pC)  is  composed  of  scalar  quantities  of  the 
constituents,  its  value  is  a  linear  combination  based  on  the  volume  frac¬ 
tions  of  fibers  and  matrix  material,  respectively, 

(pC)  -  (pC)fvf  +  (pC)mvm  (2) 

where  v^  and  vm  denote  the  fractions  of  the  fiber  volume  and  matrix 
volume,  respectively,  with  V£  +  va  -  1  naturally.  The  subscripts  f  and 
m  refer  to  the  fiber  and  matrix  constituents. 

The  effective  conductivity  k  in  Equation  (1)  is  related  to  the 
thermal  conductivities  of  the  two  different  constituents  in  terms  of  a 
parametric  relation, 

(k/k^)  -  geometry)  .  (3) 

For  heat  conduction  along  the  axial  direction  of  uni-directional  fibers. 
Equation  (3)  assumes  a  particularly  simple  form: 
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ka  “  vfkfa  +  vBkn  •  W 

For  hear  condcution  transverse  to  the  uni-directional  fibers,  the  rela¬ 
tionships  expressed  by  Equation  (3)  are  rather  involved  and  have  been 
thoroughly  analyzed  in  an  earlier  study  [3].  Equations  (1),  (2),  and 
(3)  plus  required  boundary  conditions  constitute  the  basis  for  experi¬ 
ments  from  which  temperature  data  are  to  be  analyzed  in  order  to  obtain 
the  basic  thermo-physical  property  values  of  the  constituents. 

2.  METHOD  OF  EXPERIMENTAL  DATA  ANALYSIS 

Consider  a  one-dimensional  heat  conduction  experiment  for  which 
the  physical  arrangement  is  shown  in  Figure  2.  Two  identical  test  speci¬ 
men  blocks  are  located  symmetrically  on  both  sides  of  a  centrally  loca¬ 
ted  heating  plate  assembly.  For  each  specimen  block,  one  surface  is  in 
contact  with  a  cooling  plate  which  keeps  the  specimen's  surface  tempera¬ 
ture  at  a  constant  value. 

During  a  transient  heat  conduction  experiment,  the  central  heater 
is  activated,  sending  one-half  of  the  electric  heat  input  to  the  speci¬ 
men  block.  Symmetry  of  the  physical  arrangement  assures  this  even  heat 
division.  By  monitoring  the  two  surface  temperatures,  one  on  the  heat¬ 
ing  surface  and  the  other  on  the  cooling  surface,  and  at  least  one  more 
temperature  in  an  interior  location  of  the  specimen  block,  there  are  ob¬ 
tained  three  temperature- time  histories  at  three  points,  illustrated  as 
points  A,  B,  and  C  in  Figure  2.  From  the  well-established  theory  of 
heat  conduction,  the  C  constitutes  the  boundary  conditions  in  solving 
Equation  (1) .  The  solution  for  point  B  is  uniquely  determined  if  the 
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Schematic  Arrangement  of  Teat  Apparetua  for  One-Dlaenaional 
Heat  Conduction. 
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thermal  diffusivity  value  (effectively  speaking)  a  Is  known.  Employing 
this  procedure  in  reverse,  the  recorded  temperature-time  variation  at 
point  B  therefore  determines  the  effective  diffusivity  value  for  the  com¬ 
posite  in  the  experiment  and  in  the  direction  of  the  heat  flow. 

Furthermore,  when  steady  state  heat  conduction  i3  reached 
resulting  in  a  linear  temperature  distribution  along  the  heat  path 
A-B-C,  an  effective  thermal  conductivity  can  be  determined  baaed  on  the 
steady-state  temperature  data  and  the  known  heat  flow  through  the  speci¬ 
men  block.  Thus,  in  a  single  transient  experiment,  it  is  possible  to 
determine  the  values  of  k  and  a,  from  which  an  effective  thermal  inertia 
(pC)  can  be  extracted. 

By  preparing  specimen  blocks  so  that  the  uni-directional  fibers 
are  either  aligned  with,  or  transverse  to  the  heat  flow  direction,  direc¬ 
tional  effective  thermal  conductivities  can  be  determined.  Similarly, 
by  using  specimen  blocks  with  the  matrix  material  alone,  the  basic 
thermal  conductivity  km  and  the  thermal  inertia  (pC)m  can  be  obtained. 
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EXPERIMENTAL  PROCEDURE  AND  EQUIPMENT 

1.  TEST  SPECIMENS 

A  number  of  graphite-epoxy  slabs  with  graphite  fibers  oriented 
unl-dlrectlonally  and  of  dimensions  6x6  inches  square,  0.5  inch  thick, 
were  used  In  this  experimental  study.  The  composite  specimens  were  pre¬ 
pared  by  the  Flight  Dynamics  Laboratory,  Composites  Facility,  under 
the  supervision  of  Mr.  Robert  Achard;  they  were  fabricated  by  a  standard 
process  starting  with  prepreg  tapes  which  were  supplied  by  Hercules, 

Inc.  The  manufacturer's  identifications  for  the  resin  and  fibers  were 
3501-5A  and  AS-1,  respectively;  and  according  to  the  supplier,  the  fiber 
content  in  the  prepreg  form  was  42(+3)  percent  by  weight  and  the  ply 
(tape)  thickness  was  0.0052  (+0.0004)  inches.  However,  a  direct  count 
using  a  photo-image  process  conducted  in  this  study,  as  outlined  in  Ap¬ 
pendix  B,  determined  that  the  volume  ratio  for  the  fibers  was  equal  to 
59.7  percent. 

Also  used  in  this  experimental  study  were  half-inch  thick  slabs 
of  pure  resin  (934)  supplied  by  Fibrlte  Corp.  The  resin  in  these  slabs 
was  of  essentially  the  same  chemical  composition  as  that  in  the  compo¬ 
site  slab. 
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TEST  CELL 


The  test  cell  consisted  of  three  basic  elements:  a  centrally  lo¬ 
cated  heater  assembly  and  two  cooling  plates.  Between  the  heater  assem¬ 
bly  and  the  cooling  plates,  two  identically  constructed  specimen  blocks 
were  Inserted  so  that  symmetry  of  the  geometrical  arrangement  assured 
that  half  of  the  heater  input  traveled  through  one  test  block.  A  sche¬ 
matic  diagram  has  been  presented  in  Figure  2. 

The  heater  assembly  was  composed  of  two  aluminum  surface  sheets, 
either  6x6  square  or  3  x  6  rectangle,  and  0.03  inches  thick;  the  two 
surface  sheets  encapsulated  a  foil  heater  of  the  same  dimensions  as 
those  of  the  surface  sheets.  The  heater  of  3  x  6  inches  square  was  used 
in  conjunction  with  configuration  X  for  the  heat  flow  experiments  along 
the  fiber  axis  direction.  This  configuration  and  the  other  two  are 
illustrated  in  Figure  3.  The  heater  of  6  x  6  inches  rectangle  was  used 
for  configurations  T.  The  foil  heater  was  a  thin  film-type  con¬ 
taining  electric-conducting  films  arranged  in  a  serpentine  pattern  and 
insulated  by  a  Kapton  coating.  The  overall  thickness  of  the  film  heater 
was  0.007  Inches  and  the  heaters  had  a  maximum  operating  temperature  of 
350°F.  They  were  supplied  by  Electro-Film  Co.,  Los  Angeles,  California. 

The  cooling  plates  located  at  the  two  ends  of  the  test  apparatus 
were  6x6  square,  made  of  two  aluminum  surface  sheets  enclosing  a 
series  of  zigzag  channels  through  which  cooling  water  was  passed.  The 
temperature  of  the  cooling  water  was  maintained  by  re-circulating  it  to 
and  from  an  external  tank  with  ice  as  the  heat  sink. 
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Figure  3.  Test  Specioen  Block  Orientations. 
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TEST  SPECIMEN  ORIENTATION 


For  determining  the  thermal  properties  of  the  carbon/epoxy  com* 
posites,  a  number  of  identical  specimen  slabs  were  stacked  upon  each 
other  to  form  or  simulate  a  single  test  block.  Three  different  arrange¬ 
ments  were  used  resulting  in  three  different  orientations  of  the  carbon 
fibers  relative  to  the  heat  flow  direction.  They  are  identified  in 
Figure  3  as  configurations  X,  Y,  and  Z. 

In  order  to  evaluate  the  thermal  properties  of  the  composites 
when  heat  flow  travels  along  the  fiber  axis,  configuration  X  was  used. 
For  each  of  the  two  identical  test  blocks,  six  slabs  were  used  to  form 
a  block  of  6  x  6  x  3  inches  as  shown  in  Figure  3.  The  heating  area  was 
therefore  3x6  and  faced  a  heating  plate  of  the  same  dimension. 

On  each  of  the  surfaces  of  the  heater  assembly  and  the  cooling 
plate,  three  thermocouples  were  centrally  located  to  monitor  the  surface 
temperature  variations  of  the  test  block,  corresponding  to  locations  A 
and  C  in  Figure  2.  With  test  block  surfaces  machined  flat  and  with  a 
thermal  grease  (Thermal  Paste  No.  OT-201,  Omega  Engineering,  Inc.,  Stam¬ 
ford,  Connecticut)  at  the  interfaces,  the  measured  temperatures  were 
equal  to  those  of  the  specimen  block  surfaces.  For  configurations  X  and 
Z,  heat  conduction  was  parallel  to  the  plane  of  the  slabs,  thermal 
grease  was  not  used  at  the  inter-slab  surface,  otherwise  the  heat  flow 
paths  through  the  interfacial  thermal  grease  vould  constitute  a  heat 
flow  short-circuit.  However,  for  configuration  Y,  where  the  heat  paths 
were  normal  to  the  slabs,  thermal  grease  was  used  at  the  interfaces  so 
that  the  entire  assembly  of  slabs  became  a  single  thermal  block. 


For  monitoring  the  interior  temperature-time  variations  of  the 
simulated  specimen  block  of  configuration  X,  three  thermocouples  were 
embedded  on  the  center  line  of  the  middle  slab.  They  were  located  in 
the  middle  of  the  heat  flow  path  and  at  1  inch,  respectively,  from  the 
hot  and  cold  ends.  With  x*0  and  6  denoting  the  positions  of  the  hot 
and  cold  surfaces  in  that  order,  the  interior  thermocouples  were  loca¬ 
ted  at  x*l,  3,  and  3  inches  from  the  surfaces. 

For  configuration  Y  in  Figure  3,  arranged  to  determine  the  ther¬ 
mal  conductivity  transverse  to  the  fibers,  four  specimen  slabs  (on  each 
side)  were  thermally  joined  to  form  a  single  block.  A  heater  assen&ly 
of  6  x  6  was  used,  and  interface  thermocouples  at  x  ■  0.5,  1,  1.5  in¬ 
ches  from  the  heating  surface  were  centrally  located  with  respect  to 
the  slabs  to  monitor  the  interior  temperature  variations. 

In  order  to  complete  the  three  principal  values  of  the  thermal 
conductivities,  configuration  Z,  as  illustrated  in  Figure  3,  was  used 
to  determine  the  second  transverse  effective  conductivity.  For  each 
test  block,  six  slabs  were  stacked  to  form  a  single  test  block;  and, 
as  in  configuration  X,  the  heat  flow  path  was  confined  to  the  plane  of 
the  slabs  and  no  thermal  grease  was  used  in  the  interfaces  between  slabs. 
Five  Interior  thermocouples  were  embedded  on  the  surface  of  the  middle 
slab;  they  were  at  x  ■  1,  2,  3,  4  and  5  inches  from  the  heating  surface 
and  were  lined  up,  right  in  the  middle  of  the  block  assembly. 

For  each  of  the  specimen  orientations,  two  transient  heat  flow 
tests  were  conducted,  one  with  the  heater  assembly  at  a  constant  heat 
flow  condition  and  the  other  with  the  heat  input  regulated  by  a  con¬ 
troller  so  that  the  interface  condition  simulated  a  constant  temperature 
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condition.  In  this  way,  some  degree  of  corroboration  was  achieved  in 
the  basic  thermal  data  obtained. 

Shown  in  Figure  4  is  a  partially  assemb led  test  apparatus. 

4.  TEST  IDENTIFICATION 

For  each  of  the  three  configurations  Identified  in  Figure  3  for 
the  different  fiber  orientations  relative  to  the  heat  flow  direction, 
two  transient  tests  were  performed,  one  at  a  constant  heat  flux  and  the 
other  at  a  constant  temperature.  The  test  run  identifications  follow 
the  designations  XQ,  XT,  YQ,  etc.  In  addition,  resin  slabs  were 
stacked  in  configuration  Y  of  Figure  3  and  this  test  run  was  identified 
as  test  M.  The  following  tabulation  completes  this  experimental  range. 

Test  XQ  -  Heat  Flow  Parallel  to  Fibers,  Constant  Heat  Flux 

Test  XT  -  Heat  Flow  Parallel  to  Fibers,  Constant  Surface 
Temperature 

Test  YQ  -  Heat  Flow  Transverse  to  Fibers,  Constant  Heat 
Flux  (also  transverse  to  slabs) 

Test  YT  -  Heat  Flow  Transverse  to  Fibers,  Constant  Surface 
Temperature  (also  transverse  to  slabs) 

Test  ZQ  -  Heat  Flow  Transverse  to  Fibers,  Constant  Heat 
Flux  (parallel  to  slab) 

Test  ZT  -  Heat  Flow  Transverse  to  Fibers,  Constant  Surface 


Temperature  (parallel  to  slab) 
Test  M  -  Heat  Flow  through  Resin  Slabs 


H 
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IV 


ANALYSIS  OF  TRANSIENT  TEMPERATURE  DATA 

1.  TEMPERATURE-PROFILE  MATCHING 

As  outlined  in  Section  II,  the  method  of  temperature  profile 
matching  was  used  to  analyze  and  to  extract  the  thermal  dif fusivities 
of  the  composite  along  the  three  principal  axes  of  heat  flow.  Con¬ 
ceptually,  the  method  was  based  on  the  fact  that  the  measured  tempera¬ 
ture-time  history  in  the  interior  of  the  composite  specimen  block  (as¬ 
sumed  to  be  a  complete  solid  block)  would  match  exactly  with  the  solu¬ 
tion  to  the  one-dimensional  heat  flow  equation,  provided  that  the 
experimentally  recorded  surface  temperature  variations  (at  the  heating 
surface  and  the  cooling  surface)  were  used  as  the  boundary  conditions 
to  the  differential  equation,  and  provided  that  a  correct  value  for  the 
thermal  diffuslvlty  was  used  in  the  governing  equation. 

In  order  to  determine  the  correct  effective  thermal  diffusivity 
for  a  particular  test  orientation,  i.e.,  configurations  X,  Y,  or  Z  as 
indicated  in  Figure  3,  the  reverse  to  the  above  described  matching  was 
employed.  In  other  words,  the  one-dimensional  heat  conduction  equation 
was  solved  with  the  experimentally  obtained  surface  temperature  varia¬ 
tions  as  the  boundary  condition,  but  with  thermal  diffusivity  value  a  in 
Equation  (1)  varied  so  that  for  each  diffusivity  value  assumed  in  the 
solution  scheme,  there  would  be  a  temperature  field  as  a  function  of 
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position  and  time  for  the  interior  of  the  composite  body  in  question. 
With  the  interior  temperature-time  histories  experimentally  obtained  at 
several  stations,  it  would  be  therefore  possible  to  decipher  which 
thermal  diffusivity  value  best  represents  the  thermal  characteristic  of 
the  composite.  To  simplify  the  matter,  a  single  point  in  the  interior 
was  chosen.  It  should  not  be  on  one  of  the  two  boundary  surfaces  and 
not  too  close  to  them,  for  on  either  one  of  the  two  surfaces,  the  calcu¬ 
lated  temperatures  would  be  exactly  equal  to  the  measured  values,  which 
were  used  as  the  boundary  conditions.  A  logical  choice  was,  therefore, 
a  point  half-way  between  the  heating  and  cooling  surfaces. 

In  searching  for  the  best  thermal  diffusivity  value,  so  that  the 
calculated  temperature  variation  (with  time)  at  the  midway  point  best 
described  the  experimentally  obtained  temperature  variation,  a  quantita¬ 
tive  criterion  had  to  be  established,  since  it  was  unreasonable  to  ex¬ 
pect  that  these  two  variations  would  have  fallen  on  top  of  each  other. 

A  least-squares  error  was  therefore  chosen  as  the  criterion  to  Judge  the 
best  thermal  diffusivity  value.  This  method  is  illustrated  by  consider¬ 
ing  three  tentative  values  of  the  thermal  diffusivities  a^,  02,  and  03. 
Their  values,  in  conjunction  with  the  surface  boundary  conditions,  pro¬ 
duced  three  temperature-time  history  curves  for  a  reference  point  in  the 
composite.  In  addition,  there  was  an  experimentally  recorded  tempera¬ 
ture  variation,  which  served  as  the  reference  to  judge  the  values  of  the 
thermal  diffusivities.  The  four  temperature  variations  are  illustrated 
in  Figure  5. 

According  to  the  method  of  least  squares,  an  Integrated  error  is 
defined  by  the  general  formulation: 
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where: 

t0  ”  cut-off  time,  i.e.,  the  end  of  the  transient  period 

Ta  -  temperature  variation  at  the  reference  point  from  one¬ 
dimensional  heat  conduction  equation  based  on  a. 

T  -  experimental  temperature  variation. 

exp 

Hence,  for  each  a-value  assumed,  there  was  an  associated  error  integral 
value  E  and  a  functional  relationship  can  be  established  between  them 
and  is  illustrated  in  Figure  6.  It  follows  that  the  minimum  value  of 
E  located  the  best  value  for  the  thermal  diffusivity. 

2.  ONE-DIMENSIONAL  HEAT  CONDUCTION  SOLUTION  (FINITE-DIFFERENCE) 

To  solve  the  one-dimensional  heat  conduction  Equation  (1)  with 
the  experimentally  determined  surface  temperature  variations  as  the  in¬ 
put  boundary  conditions,  it  is  entirely  possible  to  start  with  the 
exact,  analytical  solution  with  the  effective  thermal  diffusivity  a  as  a 
parameter  in  the  solution  expression.  However,  solutions  of  this  type 
are  usually  represented  by  an  infinite  series  in  conjunction  with  the 
Duhamel's  integrals,  which  are  not  easy  to  evaluate. 

Instead  of  the  exact-solution  approach,  a  finite  difference  meth¬ 
od  of  solving  the  heat  conduction  was  used.  An  explicit  time-marching 
technique  was  adopted,  with  a  stability  criterion  of: 

a(At/Ax^)  <  0.5  .  (6) 
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Figure 


Figure 


.  Illustration  of  Least-squares  Method  for  Temperature-profile 
Matching. 


.  Illustration  of  Least-squares  Error  With  Best  Diffusivity 
Value. 
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This  rendered  the  simplest  computer  routine  possible  in  order  to  obtain 
the  temperature  variations  throughout  the  composite  domain.  Eleven  to 
fifteen  space  divisions  were  used  which  were  judged  to  have  sufficient 
accuracy  in  the  time-marching  scheme. 

3.  PRELIMINARY  ESTIMATE  OF  THERMAL  DIFFUSIVITY 

In  obtaining  the  finite-difference  solutions  with  various  values 
of  the  thermal  diffusivity,  it  was  deemed  desirable  to  have  an  approxi¬ 
mate  idea  as  to  where  the  range  of  the  value  for  a  lay.  This  was  accom¬ 
plished  by  first  interpreting  the  heat  conduction  equation  in  the  fol¬ 
lowing  form: 

a  -  (3T/3t)/(32T/3x2)  (7) 

Thus,  at  a  reference  point  x  (half-way  between  the  heating  and  cooling 
surfaces) ,  the  ratio  of  the  time  derivative  to  the  space  derivative 
could  give  an  order-of-magnitude  estimate  of  the  diffusivity  value.  In 
principle,  this  method  would  give  an  exact  value;  in  practice,  the  er¬ 
rors  in  obtaining  the  derivatives  from  the  experimental  values  were  too 
large  to  have  any  credence  to  the  results. 

To  implement  this  preliminary  assessment,  consider  three  experi¬ 
mentally  obtained  temperature- time  curves  at  three  positions  xj.,  X2 ,  X3, 
equal  distance  apart  from  one  another,  as  illustrated  in  Figure  7.  At 
the  reference  position  X2,  the  time  derivative  and  space  derivative  are 
respectively  given  by  simple  formulas  and  their  ratio  gives  the  value  of 
the  thermal  diffusivity  at  u: 
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[(TJ  -  T2)/(Ti  +  T3  -  2T2)][(Ax)2/At] 


(8) 


Evaluation  of  the  right-hand  side  member  offers  some  practical  diffi¬ 
culty;  at  very  small  times,  these  four  temperatures  are  simply  too  close 
to  one  another  to  give  a  reasonable  resolution;  at  large  times,  however, 
as  the  temperature  field  tends  toward  a  steady-state  situation,  the  ra¬ 
tio  becomes  indeterminate.  Only  at  intermediate  times  does  the  above 
formula  give  a  reasonable  indication  as  to  where  a  lies  approximately. 

In  this  way,  this  method  provided  a  starting  point  to  search  for  the 
correct  diffusivity  value  by  the  method  of  temperature-profile  matching, 
which  will  be  discussed  next. 

4.  ANALYSIS  OF  TEMPERATURE  RECORD 

For  all  the  test  runs  conducted,  the  treatment  of  the  temperature¬ 
time  record  was  essentially  similar.  The  basic  steps  are  described  in 
this  section  for  test  run  XQ  which  refers  to  configuration  X,  Figure  3, 
i.e.,  heat  flow  was  along  the  fiber  axis  and  the  boundary  condition  at 
the  heating  surface  was  that  of  constant  heat  flux.  Naturally,  the 
actual  surface  temperature  variation  measured  was  used  as  the  boundary 
condition. 

By  way  of  describing  the  experimental  results  and  conditions  in 
detail,  other  experimental  runs  were  also  covered.  Cooling  water  at 
51.7°F  (+  0.6°F  variation)  was  used  in  the  cooling  plates  which  brought 
the  entire  assembly  to  an  initial  temperature  of  51.7°F  after  5  hours 
of  cooling  period.  Transient  heat  flow  commenced  by  turning  on  the 
heater  switch  so  that  a  total  of  100  watts  (+  2  variation)  heat  in¬ 
put  was  maintained.  Symmetry  of  the  test  arrangement  assured  that  a 


cons  cant  surface  flux  of  50  watts  (+  1  variation)  was  imparted  to  the 
specimen  surface.  Temperatures  recorded  were  those  at  the  heating  sur¬ 
face,  the  cooling  surface,  and  three  intermediate  points  between  the  two 
boundaries.  The  intermediate  points,  i.e.,  those  thermocouples  located 
inside  the  composite,  were  1,  3,  and  5  inches  from  the  heated  surface. 
Figure  8  illustrates  the  five  temperature- time  histories  obtained  from 
this  test  run. 

The  data  displayed  in  Figure  8  were  treated  as  follows:  first, 
a  crude  estimate  of  the  thermal  dlffusivlty  was  made  based  on  the  formu¬ 
lation  in  Section  4.3.  At  a  selected  time  instant,  the  three  tempera¬ 
tures  were  those  at  positions  of  x  ■  1,  3,  and  5  inches,  respectively, 
i.e..  Ax  »  (1/6)  ft.;  the  temperature-time  gradient  was  obtained  by 
subtracting  the  two  successive  measurements  of  the  temperature  at  x  ■  3 
Inches.  In  this  manner,  an  a-value  was  obtained.  By  repeating  the  pro¬ 
cedure  at  different  time  Instants,  a  record  of  "variation"  of  thermal 
dlffusivlty  with  time  was  obtained;  this  particular  run,  XQ,  is  shown 
in  Figure  9.  The  "variation"  of  the  dlffusivlty  values  were  of  course 
caused  by  the  inaccuracy  inherent  in  the  formulation  of  Equation  (8). 
However,  the  record  gave  an  Indication  as  to  the  limits  in  which  the 

correct  value  lay;  it  was  within  a  lower  and  upper  limit  of  0.25  x  10“^ 

—  5  2 

and  1.5  x  10  ft  /sec.  From  these  limiting  values,  the  finite-differ¬ 
ence  equation  was  invoked  to  determine  the  error  integral  E,  defined  by 
Eq tuition  (5),  with  the  experimentally  determined  temperatures  at  x  ■  3 
inches  as  the  reference  temperatures.  By  varying  the  dlffusivlty  value, 
various  error  Integrals  were  obtained,  and  the  minimum  value  of  the  lat¬ 
ter  was  found  at  a  dlffusivlty  value  of  a  *  3.3  x  10"^  ft^/sec.  Figure  10 
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(Heat  Flow  Parallel  to  Fibers). 
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3.6 


Figure  10.  Least  Squares  Error  Distribution  for  Case  XQ 


shows  Che  computed  distribution  of  the  error  integral  E  at  various  as¬ 
sumed  values  of  a.  For  test  run  XQ,  a  value  of  3.3  x  10“^  for  effective 
thermal  diffusivlty  in  the  direction  along  the  fiber  axis  was  obtained. 

At  the  conclusion  of  the  transient  period,  a  steady-state  tempera¬ 
ture  distribution  was  obtained,  which  was  used  to  calculate  the  effec¬ 
tive  thermal  conductivity  along  the  fiber-axis  direction.  Such  a 
steady-state  temperature  record  was  indicated  in  Figure  11,  which  also 
showed  the  temperatures  at  selected  time  instants.  The  effective  ther¬ 
mal  conductivity  turned  out  to  be  3.63  BTU/hr-°F-f t. 

Having  obtained  the  effective  diffusivlty  value,  one  would  like 
to  compare  it  with  the  calculated  temperatures  from  the  one-dimensional 
differential  Equation  (1),  using  the  effective  diffusivlty  of  3.3  x  10”^ 
ft  /sec.  For  test  run  XQ,  three  temperature  records  at  x  ■  1,  3,  and  5 
inches,  respectively,  were  available;  their  temperature  variations,  to¬ 
gether  with  the  calculated  variations,  are  shown  in  Figure  12.  The  cal¬ 
culated  values  at  x  -  1  inch  were  above  the  experimental  while  those  at 
x  “  5  inches  were  below  the  measured  data.  For  the  reference  position 
of  x  ■  3  Inches,  agreement  appears  satisfactory. 

The  individual  diffusivlty  and  conductivity  values  from  all  test 
runs  were  extracted  in  a  similar  fashion  and  were  summarized  in  Table  1. 

A  noteworthy  feature  of  the  data  was  the  fact  that  the  effective  thermal 
conductivity  for  heat  flow  parallel  to  the  fiber  axis  was  one  order  of 
magnitude  above  the  corresponding  value  for  heat  flow  transverse  to  the 
fiber  axis. 

The  transient  temperature  records  for  all  test  runs  as  well  as 
comparisons  of  the  predicted  values  at  the  interior  points  based  on  the 
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Recorded  Temperatures 


Figure  11.  Selected  Temperature-time  Records  for  Case 


Calculated  Temperature 


arlsons  of  Calculated  Temperature  Responses  Versus  Experimental 
ies.  Case  XQ. 


TABLE  1 


SUMMARY  OF  EFFECTIVE  CONDUCTIVITIES  AND  EFFECTIVE  DIFFUSIVITIES 


Test  Run 

Effective  Conductivity 
k  (BTU/hr-°F-ft) 

Effective  Diffusivity 
a  (ft2/sec) 

Remarks 

XQ 

XT 

3.59 

3.66 

3.3  x  10*  5 

3.0  x  10'5 

Heat  flow 
along  fiber 
axis. 

Configuration 

X,  Figure  3. 

Average 

3.63 

3.15  x  10"5 

YQ 

0.452 

2.0  x  10“J 

Heat  flow 

0.459 

2.7  x  10"6 

transverse  to 

fibers  and 

soeclmen  slabs 

YT 

0.431 

2.2  x  10"° 

Configuration 

0.451 

2.8  x  10“6 

Y,  Figure  3 

Average 

0.448 

2.4  x  10*6 

ZQ 

0.700 

3.19  x  10“6 

Heat  flow 

transverse  to 

fibers,  paral- 

lei  to  slabs 

ZT 

0.590 

2.71  x  10‘ft 

Configuration 

Z,  Figure  3 

Average 

0.645 

2.95  x  10~6 

M 

0.136 

8.5  x  10~7 

Heat  flow  nor- 

mal  to  resin 
slabs ,  con¬ 
stant  heat 
flux 

extracted  thermal  diffusivity  values  are  contained  in  Aopendix  A,  which 
also  Includes  the  temperature- time  records  for  test  run  XQ  illustrated  in 
this  section. 

Some  comments  are  in  order  with  regard  to  the  effective  thermal 
conductivities  in  the  directions  parallel  and  transverse  to  the  fiber 
axis:  the  values  reported  are  for  an  average  temperature  of  150°F  and  the 
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axial  value  is  8.10  times  as  large  as  the  lower  one  of  the  two  trans¬ 
verse  values.  In  Reference  [4],  effective  conductivities  for  SO  percent 
uni-directional  fibrous  composites  are  listed  as  7.05  and  0.55  BTU/ 
hr-°F-ft  (at  a  temperature  of  350°K)  for  the  axial  and  transverse  values, 
respectively.  Their  ratio  is  13.  A  difference  of  the  volume  ratio  0.50 
from  the  value  used  in  this  study  may  be  partially  responsible.  However, 
the  effective  value  of  7.05  reported  in  [4]  is  at  considerable  variance 
with  the  value  of  3.63  given  in  Table  1;  the  difference  may  be  attributed 
to  the  different  tow  structures  and  the  different  carbonization  tempera¬ 
tures  used  in  their  respective  processes. 
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V 


THERMO-PHYSICAL  PROPERTIES  OF  THE  CONSTITUENTS 

Having  obtained  the  directional  values  for  the  effective  thermal 
conductivities  and  diffusivities ,  the  next  step  was  to  deduce  basic 
thermo-physical  properties  of  the  constituent  materials.  The  scalar  pro¬ 
duct  (pC)  for  the  carbon/epoxy  composite  can  be  calculated  from  data 
obtained  for  each  of  the  three  configurations  X,  Y,  and  Z.  Starting 
with  the  definition  pC  *  k/a,  the  effective  thermal  inertia  values  were 
obtained  for  each  configuration;  these  values  should  conceptually  agree, 
but  have  some  variation  as  noted  below: 


From  Configuration  X: 

(PC) 

-32.0  BTU/ft3-°F 

Y: 

(PC) 

-  51.85  BTU/ft3-°F 

Z: 

(PC) 

-  60.73  BTU/ft3-°F 

Average  48.19 

From  Configuration  M:  (pC)m  -  44.4  BTU/ft3-°F 
(pure  resin  slabs) 


It  should  be  pointed  out  that  the  values  for  thermal  inertia  (oC)  for  the 
composite  specimens  should  agree  with  engineering  accuracy  based  on  the 
data  from  all  three  configurations  X,  Y,  and  Z.  However,  the  data  based 
on  configuration  X  revealed  some  Inherent  inaccuracy  in  the  reduction  of 
the  data  to  obtain  the  effective  conductivity  and  effective  diffuslvity. 
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For  Chis  reason,  less  credence  was  placed  on  the  (pC)-value  from  that 
configuration.  It  must  also  be  recognized  that  the  product  term  (pC) 
was  obtained  by  dividing  k  by  a,  both  of  which  were  deduced  from  the  ex¬ 
perimental  data.  In  obtaining  the  ratio,  errors  in  the  individual  nu¬ 
merical  values  were  amplified.  Taking  the  average  of  the  three  (pC)- 
values,  an  accepted  result  was: 

(pC)  ■  48.2  BTU/ft3-°F,  for  composites. 

Further  deductions  for  the  basic  thermo-physical  properties  of  the  con¬ 
stituents  are  discussed  below. 

1.  DENSITY 

Manufacturer's  specifications  for  the  prepreg  that  went  into 
graphite /epoxy  specimen  gave  the  density  values  of  0.0660  and  0.0460 
lb^/in  for  fiber  and  resin,  respectively.  Graphite-epoxy  specimen 
slabs  and  resin  slabs  were  weighed  and  their  volumes  noted  in  connection 
with  the  present  study;  from  these  observations,  there  resulted  the  fol¬ 
lowing  density  values: 

For  composite  (0.597  fiber,  0.403  resin  by  volume) 
p  -  0.0578  lb  /in3 

m 

For  resin 

p  ■  0.0462  lb_/in3  . 

B  m 

Since  the  density  of  a  composite  is  related  to  the  constituent  densities 
by  the  relation: 
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P  "  pfvf  +  pmvm  ’  (9) 

its  fiber  density  can  be  computed  from  Equation  (9),  with  Vj  “  0.597, 
and  v_  ■  0.403.  This  resulted  in: 

pj  »  0.0656  lbm/in3  . 

2.  SPECIFIC  HEAT 

For  the  resin  matrix,  the  values  of  (pC)  and  p  gave  the  specific 

IS  IS 

heat  for  resin  as: 

Cm  -  0.556  BTU/lbm-°F  . 

For  fiber  specific  heat,  the  calculations  were  not  as  direct  as 
for  the  resin  matrix.  First,  the  thermal  inertia  term  (pC)  for  the  com¬ 
posite  was  related  to  the  inertia  terms  of  the  constituents  by  the  rela¬ 
tion: 


(pC)  -  (pC)fvf  +  (pC)Bvm  .  (10) 

The  best  values  for  (pC)  of  the  composite  and  for  (pC)m  of  the  resin  ma¬ 
trix  were  inserted  into  Equation  (10),  with  v^  and  vB  at  0.597  and  0.403, 
respectively.  The  product  (pC)^  for  the  fiber  was  determined  to  be: 

(pC)f  -  50.8  BTU/ft3-°F  . 
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Since  p£  has  been  determined  as  0.065b  lb/in^,  Che  specific  heat  for 
carbon  fibers  was  therefore: 

C£  -  0.448  BTU/lbm  -  °F  . 

3.  THERMAL  CONDUCTIVITY 

For  the  resin  matrix,  which  was  assumed  isotropic,  the  thermal 
conductivity  was  listed  in  Table  1  as  0.136  BTU/(hr-°F-f t) . 

For  the  fibers  in  the  composite,  there  were  three  principal  va¬ 
lues.  From  the  data  obtained  with  heat  flow  parallel  to  the  uni¬ 
directional  fibers,  the  fiber  conductivity  along  the  fiber  axis  can  be 
deduced  based  on  a  parallel  heat  flow  circuit  expressed  by  Equation  (4). 
From  Table  1,  the  effective  conductivity  for  heat  flow  along  the  uni¬ 
directional  fiber  is  given  as  3.63  BTU/(hr-°F-ft) ;  since  has  been  de¬ 
duced,  therefore  the  axial  conductivity  for  fibers  can  be  obtained. 

Icfa  -  5.99  BTU/(hr-°F-ft) 

along  the  fiber  axial  direction. 

For  the  transverse  conductivity  across  the  fiber  axis,  the  method 
of  deduction  was  not  as  simple. 

There  were  two  measured  values  for  the  effective  conductivities 
when  heat  flow  is  transverse  to  the  uni-directional  fibers.  Normal  to 
the  slab-specimen  plane,  the  effective  value  obtained  was  0.448  BTU/ 
hr-°F-ft;  in  the  slab-specimen  plane  the  effective  value  is  0.645.  Inso¬ 
far  as  the  dispersion  pattern  shown  in  Figure  1  exhibits  no  preferential 


direction,  Che  difference  between  these  two  values  were  the  results  of 
heat  leakage  to  Che  side.  In  examining  the  experimental  arrangements, 
the  data  obtained  with  heat  flow  normal  to  the  slab-plane  was  considered 
more  reliable  since  thermal  grease  was  used  to  minimize  the  contact  re¬ 
sistances  at  inter-slab  locations  among  other  factors  considered. 

In  Che  analysis  of  [3],  the  hexagonal  dispersion  pattern  vas 
shown  to  produce  equal  effective  conductivities  regardless  of  the  heat 
flow  direction,  therefore  the  question  is  now:  what  is  the  ratio  of  the 
fiber  (transverse)  conductivity  to  the  matrix  (isotropic)  conductivity 
for  a  fiber-volume  ratio  of  0.597  such  that  the  ef fective-to-matrix  con¬ 
ductivity  ratio  is  3.29  (*■  0.448/0.136)?  By  using  a  built-in  subroutine 
(STAG),  the  ratio  turned  out  to  be  17.8  or  the  transverse  conductivity  of 
the  fibers  is  17.8  x  .136  ■  2.42  BTU/hr-°F-ft2/ft.  Hence,  the  fibers 
were  shown  in  the  above  analysis  to  have  two  conductivities  as  follows: 

(i)  Along  the  fiber  axis:  Kfa  ■  5.99  BTU/hr-°F-f t2/f t 
(ii)  Transverse  to  the  fiber  axis:  ■  2.42  BTU/hr-°F-f t2/f t 

There  were  very  few  data  reported  in  the  literature;  the  only  pub¬ 
lication  which  listed  the  axial  conductivity  for  carbon  fibers  was  that 
of  Larsen  [5],  who  reported  a  value  of  2.5  3TU/hr-°F-f t2/f t  which  was 
about  one-third  of  the  value  measured  in  this  study.  It  is  of  interest 
to  note  that  the  isotropic  conductivity  of  resin  matrix  was  also  given  by 
Larsen  in  Reference  [5],  as  0.133  BTU/hr-°F-ft2/ft  which  was  in  excellent 
agreement  with  the  value  of  0.136  measured  in  this  experimental  work. 
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Another  value  for  the  axial  conductivity  of  carbon  fibers  was  de¬ 
duced  from  the  data  found  in  Reference  [4].  An  axial  effective  conducti¬ 
vity  for  a  50  percent  carbon-epoxy  composite  (uni-directional)  was  listed 
as  0.117  tf/Cm-°K,  or  6.76  BTU/(hr-°F-f t)  for  a  temperature  of  150°F. 

Using  Equation  (4)  to  synthesize  for  the  fiber  axial  conductivity,  a  va¬ 
lue  of  k^a  ■  13.4  BTU/(hr-°F-ft)  was  obtained.  (In  doing  so,  kj,  was  set 
at  0.136,  which,  even  if  inaccurate,  hardly  affected  the  result.)  From 
these  sources,  then,  the  axial  conductivity  for  carbon  fibers  ranged  from 
2.5  BTU/(hr-°F-f t)  found  in  Reference  [5]  to  13.4  BTU/(hr-°F-ft)  based  on 
the  data  of  Reference  [4],  while  the  value  determined  in  this  investiga¬ 
tion  turned  out  to  be  6.0  BTU/(hr-°F-f t) .  The  variation  was  considerable 
indeed. 
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CONCLUSIONS 


The  purposes  of  this  experimental  investigation  were  twofold:  one 
was  to  explore  the  thermal  physical  properties  of  a  carbon/epoxy  compo¬ 
site,  a  frequently  used  material  in  aerospace  structures,  and  the  other 
was  to  develop  a  reliable  measurement  method  to  determine  these  proper¬ 
ties  by  using  carbon/epoxy  materials  as  an  example.  By  employing  a 
transient  heat  flow  technique,  the  effective  thermal  conductivities  and 
diffusivities  along  three  principal  directions  were  determined  for  car¬ 
bon/epoxy  composites  with  uni-directional  fibers;  uni-directional  orient¬ 
ation  was  chosen  because  it  afforded  a  simpler  geometry  of  the  dispersion 
pattern  of  the  fibers  than  other  dispersion  prescriptions. 

The  thermal  conductivities  of  the  constituents  -  the  fibers  and 
the  matrix  -  of  the  composites  were  determined  by  first  establishing  the 
dispersion  pattern  of  the  fibers  through  microscopic  photographic  examin¬ 
ations  of  the  specimen  cross-section,  illustrated  in  Figure  1;  and  by 
using  the  relation  between  the  ef fective-to-matrix  conductivity  ratio 
and  the  fiber-tcr-matrix  conductivity  ratio  obtained  in  Reference  [l]  or 
AFWAL-TR-80-3012. 

With  regard  to  the  micro-structure  of  the  carbon/epoxy  composites 
prepared  from  plies,  microscopic  examinations  indicated,  generally  speak¬ 
ing,  a  uniform  dispersion  pattern  very  much  resembling  a  staggered 


arrangement  into  hexagonal  arrays.  An  unusual  feature  was  the  excess 
matrix  (resin)  materials  which  congregated  near  the  interfaces  between 
plies  and  formed  a  canal-like  border  of  one-tenth  the  thickness  of  an 
entire  ply.  Thus,  in  the  cross  section  plane,  the  structural  pattern 
was  not  homogeneous  macroscopically  and  this  factor  contributed  to  the 
difference  between  the  two  transverse  effective  thermal  conductivities, 
and  led  undoubtedly  to  different  elastic  properties  in  these  two  trans¬ 
verse  directions  as  well.  Further  investigations  are  needed  to  clarify 
the  observed  discrepancy. 

The  method  of  synthesizing  the  transverse  effective  thermal  con¬ 
ductivities  to  obtain  fiber-to-matrlx  conductivity  ratios  was  based  on 
an  exact  mathematical  analysis  [3];  in  the  range  of  this  experimental  in¬ 
vestigation — fiber/epoxy  composites  with  0.60  fiber  content  by  volume — 
the  relationship  was  such  that  a  large  change  in  the  fiber/epoxy  thermal 
conductivity  ratio  produced  a  relatively  small  change  in  the  effective 
conductivity  ratio.  Since  the  Inverse  process  was  used  for  synthesis, 
i.e.,  from  an  effective  value  to  the  constituent  value,  the  process  was 
therefore  somewhat  insensitive,  and  for  further  work,  composites  with  low 
volume  ratios  for  fibers  are  more  desirable  from  the  sensitivity  stand¬ 
point. 

The  measurements  conducted  in  this  study  were  limited  by  instru¬ 
mentation  problems  to  a  level  of  250°F.  Since  fiber/epoxy  composites 
are  known  to  exhibit  temperature-dependence  of  their  thermal  properties, 
future  efforts  are  required  to  extend  the  temperatures  to  a  more  practi¬ 
cal  range. 
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APPENDIX  A 


TEMPERATURE  DATA  AND  COMPARISONS  WITH  PREDICTIONS 


45 


In  the  following  pages,  the  temperature- time  records  are  present¬ 
ed  for  the  test  runs  on  composite  specimens  in  configurations  X,  Y,  and 
Z  as  well  as  on  matrix  (resin  only)  specimens.  For  clarification  of 
these  configurations,  reference  is  made  to  Figure  3,  page  IS.  Two  tests 
were  performed  for  each  configuration,  one  with  the  heating  surface  at  a 
constant  heat  flux  condition  and  the  other,  at  a  constant  surface  temper¬ 
ature  condition. 

Experimental  data  shown  in  dashed  lines  are  compared  with  theore¬ 
tical  predicted  values  which  are  based  on  the  diffusivity  values  obtained 
by  a  best  fit  of  the  experimental  curve  with  the  theoretical  prediction 
for  the  central  point,  midway  between  the  heating  and  cooling  surface. 
These  two  measured  surface  temperatures  were  used  as  the  boundary  condi¬ 
tions  in  obtaining  the  interior  temperature- time  histories.  The  best-fit 
thermal  diffusivity  values  for  each  test  run  are  listed  in  Table  1,  page 
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gure  A-l.  Temperature  Profile  Matching  Results,  Configuration  X,  Constant  Heat 
Flux  (XQ). 


ature 


Figure  A- 2. 


Tenpereture  Profile  Matching  Results 
etant  Surface  Tenperature  (XT). 


Configuration  X,  Con- 


Figure  A-4.  Temperature  Profile  Matching  Results,  Configuration  Y,  Con¬ 
stant  Surface  Temperature  (YT1) . 
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Experimental 


Figure  A-5,  Temperature  Profile  Matching  Results,  Configuration  Y,  Constant  Heat 
Flux  (YQ2) . 


Figure  A-6.  Temperature  Profile  Matching  Results,  Configuration  Y,  Con¬ 
stant  Surface  Temperature  (YT2). 
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Temperature, 


Figure  A— 8.  Temperature  Profile  Matching  Results,  Configuration  Z,  Con¬ 
stant  Surface  Temperature  (ZT). 
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Temperature, 


Time,  1 03  sec. 


Figure  A-10.  Temperature  Profile  Matching  Results,  Resin  Specimen,  Con¬ 
stant  Surface  Temperature  (MT). 


The  volume  ratio  of  the  graphite-epoxy  composite  used  In  this 
study  was  calculated  by  visually  counting  and  measuring  the  fibers  (in 
reality,  tows  of  fibers)  in  a  given  cross-sectional  area,  similar  to 
Figure  1.  An  enlarged  photograph  of  the  composite,  250X,  was  sectioned 
into  fifteen  equal  units  of  approximately  two-hundred  fibers  each,  and 
the  average  number  of  fibers  per  unit  was  determined.  Then  one  section 
was  enlarged,  overall  magnification  2100X,  and  one-hundred  fiber  diame¬ 
ters,  randomly  chosen,  were  measured.  The  diameter  measurement  was  re¬ 
peated  in  order  to  assuage  individual  observational  bias.  Next,  the 
diameters  (actually  the  squares  of  the  diameters)  were  given  a  statisti¬ 
cal  distribution  test,  the  x~square  test,  which  showed  that  the  quantity 
D  did  follow  a  normal  (Gaussian)  distribution.  An  average  image  diame¬ 
ter,  D  ■  0.630  inches,  was  obtained  from  the  second  enlargement  (which 
corresponded  to  an  actual  average  diameter  of  0.0030  inches*),  and  was 
therefore  assumed  to  apply  to  the  other  fourteen  counting  units. 

V  -  Nf(wD2/4A) 

*Hercule's  specifications  say  that  the  ply  thickness  is  0.005  inches,  so 
the  value  for  the  fiber  diameter  is  compatible. 
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V 


volume  ratio 


■  average  number  of  fibers  in  the  area  counted 

—2 

D  -  square  of  the  average  diameter 
A  -  counting  area 

The  entire  process  gave  a  fiber  volume  ratio  of  59.7  percent  with  a  4 
percent  standard  deviation,  i.e.,  the  volume  ratio  lies  within  the  limits 
of  55.7  to  63.7  percent  with  59.7  as  the  most  probable  range.  The  popu- 
lation  distribution  for  the  quantity  D  (D  ■  image  diameter,  inches)  for 
the  first  counting  unit  is  shown  in  Figure  B-l. 
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Figure  B— 1.  Fiber  Diameter  Population  Distribution. 


The  purpose  of  Che  experimental  apparatus  was  to  generate  and  re¬ 
cord  transient  temperature  data  present  in  the  graphite-epoxy  block.  In 
order  to  accomplish  this,  equipment  was  obtained  which  controlled  the 
temperature  at  the  heating  face  and  cooling  face,  monitored  and  recorded 
transient  temperature  data,  and  insulated  the  graphite-epoxy  test  speci¬ 
men. 

Enclosed  herein  is  the  following:  a  summary  of  the  equipment  which 
was  used  to  conduct  the  two  experiments,  a  detailed  description  of  equip¬ 
ment  items  and  experimental  procedure;  a  section  outlining  the  function 
of  the  three  equipment  groups,  i.e.,  temperature  control  system,  the  in¬ 
sulation,  temperature  measurement  system,  and  finally  a  specification 
list  of  the  individual  equipment. 

Two  boundary  conditions  were  imposed  at  the  heating  face  of  the 
graphite-epoxy  composite:  one  was  a  constant  heat  flux,  Q,  and  the  other 
was  a  constant  temperature  step,  T.  Table  C-l  summarizes  equipment  used 
in  each  experiment.  The  two  experiments  were  then  conducted  as  follows: 

Constant  Surface  Temperature  Experiment  (T) 

Equipment  items  used  for  this  part  of  the  experiment  and  their 
Interconnections  are  illustrated  in  Figure  C-l.  Temperature  signals  from 
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TABLE  C-l 


SUMMARY  OF  EXPERIMENTAL  EQUIPMENT 


Equipment 

Experiment 

Case 

Variable  Transformer 

Q 

All 

Temperature  Controller 

T 

All 

Clock 

T 

All 

Digital  Voltmeter 

Q,T 

All 

Watt  Transducer 

Q»T 

All 

6"  x  6"  Heating  Plate 

Q,T 

Y 

3"  x  6"  Heating  Plate 

Q»T 

i  m 

Thermal  Feedback  Sensor 

T 

the  test  cell  were  recorded  by  a  Fluke  Data  Logger,  and  simultaneously, 
the  cooling  plate  temperature  was  monitored  by  a  thermocouple  reader. 

For  maintaining  the  surface  temperature  at  a  pre-set  value,  a  temperature 
controller  (Chromalox  3804)  was  used  to  sense  the  surface  temperature  of 
the  specimen  block  and  adjust  the  power-supplier  accordingly.  A  timing 
device,  built  into  the  circuit,  monitored  the  duration  in  which  power  was 
supplied  to  the  heater;  in  this  way  the  Integrated  power  input  could  be 
obtained,  which  was  later  needed  for  steady-state  calculations  of  the 
thermal  conductivity  values.  The  power  supply  rating  was  displayed 
through  a  Watt  transducer  (PT  2S01F,  Bell).  The  assembled  items  are 
pictorlally  presented  in  Figure  C-2;  this  also  shows  the  cooling  tank  for 


Temperature 

Controller 


Thermocouple 

Reader 


Figure  C-2.  Assembly  of  Experimental  Equipment. 


maintaining  the  cooling  surface  temperature. 

The  sequential  steps  in  conducting  this  experiment  consisted  of 

the  following: 

1.  The  test  block  was  placed  inside  the  test  chamber,  as  illustrated 
in  Figure  4,  and  all  electrical  leads  were  connected  to  their  pro¬ 
per  terminals. 

2.  The  cooling  tank  valve  was  opened  to  cool  the  specimen  to  its 
Initial  temperature,  approximately  50°F.  This  temperature  was  as¬ 
sumed  to  have  been  reached  after  Che  temperature  gradient  in  the 
test  cell  was  less  than  or  equal  to  0.5°F. 

3.  The  temperature  control  was  then  set  so  that  the  temperature  at 
the  heating  face  was  equal  to  the  initial  temperature  of  the  test 
cell  plus  200°F.  Upon  reaching  its  initial  condition,  a  constant 
step  change  in  temperature  was  imposed  via  the  temperature  con¬ 
troller. 

4.  Temperatures  inside  the  test  cell  were  recorded  at  the  heating 
face,  and  at  three  locations  between  the  two  boundary  plates. 

3.  The  temperature  at  the  cooling  plate  was  monitored  continuously 

and  kept  within  Q.5°F  limit.  The  heating  plate  was  controlled  by 
the  temperature  controller;  its  temperature  oscillated  no  more 
than  3°F  as  a  consequence. 

6.  Data  were  recorded  at  approximately  5  to  10°F  intervals  by  the 
Data  Logger  such  that  a  smooth  temperature  profile  could  be  con¬ 
structed. 

7.  Steady  state  and  the  conclusion  of  the  experiment  were  assumed 
after  the  temperature  change  at  all  stations  was  less  than  0.2°F/ 


1000  seconds 
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8.  At  the  conclusion  of  the  experiment,  the  following  information  was 

recorded: 

(i)  initial  conditions 

(ii)  temperature  profile  at  steady  state 

(iii)  average  wattage  supplied,  Q 

(iv)  operational  time,  "on-time"  of  the  Chromolox  temperature 
controller,  tc 

(v)  total  elapsed  time  of  the  experiment  t 

o 

(vi)  fiber  orientation 

(vii)  power  required  at  steady  state,  q  ■  Q(tc/tQ) 

Constant  Flux  Experiment  (Q) 

The  equipment  items  used  in  this  experiment  were  essentially  simi¬ 
lar  to  those  employed  in  the  constant  surface  temperature  part.  The 
chief  difference  was  that  a  variable  transformer  was  used  in  the  constant 
heat-flux  experiment  in  lieu  of  the  PID-controller,  and  also,  the  timing 
device  was  disconnected.  The  apparatus  schematic  is  illustrated  in  Fi¬ 
gure  C-3.  Temperature  signals  from  the  test  cell  were  recorded  by  the 
Fluke  Data  Logger,  and  the  cooling  plate  temperature  was  monitored  simul¬ 
taneously  by  a  thermocouple  reader.  Finally,  the  amount  of  power  sup¬ 
plied  to  the  heating  plate  was  also  monitored  using  a  Watt  transducer. 

In  order  to  effect  a  200°F  temperature  gradient  between  the  heat¬ 
ing  and  cooling  plate,  the  power  setting  on  the  variable  transformer  had 
to  be  known.  This  value  of  power  setting  was  obtained  from  the  constant 
surface-temperature  experiment.  The  experiment  was  then  conducted  as  fol¬ 
lows: 
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Variable 

Transformer 


1.  The  test  block  was  assembled  Into  the  test  cell,  and  connected  to 
the  proper  equipment. 

2.  The  cooling  tank  valve  was  opened  to  cool  the  specimen  to  its  ini¬ 
tial  temperature,  approximately  50°F;  this  took  from  five  to  ten 
hours.  A  uniform,  initial  temperature  was  assumed  when  the  tem¬ 
perature  variation  in  the  test  cell  was  less  than  or  equal  to 
0.5°F. 

3.  The  variable  transformer  was  then  set  so  that  a  temperature  gradi¬ 
ent  of  200°F  would  result  across  the  graphite-epoxy  block  at  the 
conclusion  of  the  experiment. 

4.  The  temperatures  inside  the  test  cell  were  recorded  at  the  heating 
face,  cooling  face,  and  three  locations  inside  the  graphite-epoxy 
test  block. 

5.  The  temperature  at  the  cooling  plate  was  monitored  continuously 
and  kept  within  0.5°F.  The  heating  plate  was  controlled  by  the 
variable  transformer. 

6.  Data  were  recorded  at  approximately  5  to  10°F  Intervals  by  the  Data 
Logger  such  that  a  smooth  temperature  profile  could  be  constructed. 

7.  Steady  state  and  the  conclusion  of  the  experiment  was  assumed  when 
the  temperature  rise  at  all  stations  was  less  than  0.2°F/1000  se¬ 
conds  . 

8.  At  the  conclusion  of  the  experiment,  the  following  information  was 
recorded: 

(i)  initial  condition 

(ii)  temperature  profile  at  steady  state 
(ill)  average  wattage  supplied 
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(iv)  total  elapse  time  of  the  experiment 

(v)  fiber  orientation 

Temperature  Control  System 

(i)  Heating  System 

The  object  of  the  heating  system  was  to  impose  a  boundary  condi¬ 
tion  at  the  heating  plane  surface.  It  consisted  of  a  power  supply  and 
an  electrically  resistant  heating  plate.  Two  different  power  supply 
units  were  used  in  the  experiments:  one  was  a  variable  transformer  and 
the  other,  a  proportional-integral-differential  (PID)  temperature  con¬ 
troller.  The  transformer  was  used  in  those  experiments  requiring  a  con¬ 
stant  heat  flux,  Q;  and  the  PID-controller  was  used  in  the  alternate  ex¬ 
periments  for  a  temperature  step  change,  T. 

Two  heating  plates,  a  3  x  6  x  0.0625  inch  and  a  6  x  6  x  0.0625 
inch,  were  built.  The  rectangular  heating  plate,  3x6,  was  constructed 
by  sandwiching  three  1x6  inch  heating  strips  between  two  3x6  inch 
aluminum  plates.  The  strip  heaters  were  then  wired  in  parallel,  and  the 
heating  plate  was  instrumented  with  three  iron- cons tan tan  thermocouples — 
one  per  heating  strip.  The  center  thermocouple  was  used  to  monitor  the 
temperature  of  the  heating  plate  surface,  and  the  outside  thermocouples 
were  connected  in  parallel  and  used  either  as  a  feedback  signal  to  the 
PID-controller  or  as  a  check  of  the  center,  recorded  temperature.  The 
3x6  heating  plate  was  used  with  the  graphite-epoxy  composite  configura¬ 
tions  X  and  Z. 

The  6x6  heating  plate  was  built  in  the  same  manner  as  the  3x6 
plate.  A  6  x  6  electrically  resistant  heating  strip  was  placed  between 
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two  6x6  inch  aluminum  plates  and  instrumented  with  three  iron- cons tan- 
tan  thermocouples.  Two  thermocouples  were  wired  in  parallel  and  were 
used  to  record  the  temperatures  of  the  heating  plate  surface.  The  third 
thermocouple  was  used  either  as  a  signal  feedback  to  the  PID-controller 
or  as  an  extra  recorded  temperature  data  point.  The  6x6  heating  plate 
was  used  during  those  experiments  for  configuration  Y . 

(ii)  Cooling  System 

The  purpose  of  the  cooling  system  was  to  reduce  the  temperature 
of  the  graphite-epoxy  test  specimens  to  an  initial  temperature  and  to 
maintain  the  boundary  condition  at  the  cooling  plane.  The  cooling  sys¬ 
tem  consisted  of  the  following  items: 

1)  cooling  tank,  30  gallons 

2)  cooling  plates,  6x6x1  inch 

3)  catch  tank,  10  gallons 

A)  return  pump,  20  gallons  per  minute 
S)  thermocouple  indicator,  iron-cons  tan tan 

Water  was  cooled  with  ice  in  the  cooling  tank  to  approximately  50°F  and 
then  circulated  through  the  cooling  plates  into  the  catch  tank  and  back 
to  the  cooling  tank.  Constant  temperature  at  the  cooling  plate  was  main¬ 
tained  by  operator  observation.  Three  thermocouples  were  attched  to  the 
cooling  plates;  two  were  wired  in  parallel  and  connected  to  the  Data  Log¬ 
ger,  and  the  third  was  attached  to  the  thermocouple  indicator  and  moni¬ 
tored  continuously.  A  change  in  temperature  plate  of  more  than  0.5°F  was 
corrected  by  the  operator  manually. 
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Insulation 

The  Insulation  provided  an  adiabatic  condition  on  the  peripheral 
surfaces  of  the  graphite-epoxy  specimens ,  and  this  shell,  or  test  cell, 
was  then  placed  into  a  fiberglass-insulation  chamber.  The  insulation 
chamber  was  built  using  approximately  fifty,  4'  x  4'  x  2"  fiberglass 
panels.  The  chamber  itself  was  approximately  4-foot  cube,  and  the  test 
cell  was  placed  into  the  center  of  the  cube.  The  R-value  of  the  insula¬ 
tion  at  the  center  was  approximately  equal  to  90,  i.e.,  k  -  1.12  x  10-i 
BTU/(ft-hr-°F). 

Temperature  Recording  System 

For  recording  the  boundary  temperatures  and  the  transient  tempera¬ 
tures  inside  the  graphite-epoxy  composite,  the  system  consisted  of  a 
Fluke  Data  Logger  and  multi- junction,  iron-cons tan tan  thermocouples. 

Three  thermocouples  were  wired  in  parallel,  thus  forming  a  single  group 
of  parallel- junction  thermocouples.  Three  groups  of  these  were  made. 

The  three  parallel- junction  thermocouples  were  then  attached  to  one  of 
the  graphite-epoxy  specimens.  In  the  graphite-epoxy  composite,  shallow 
grooves,  approximately  1/8  x  1/8  inch  were  milled  at  1,  3,  and  3  inches 
from  the  heating  face,  and  the  three  parallel-junction  thermocouples  were 
cemented  below  the  plane  surface  near  the  geometrical  center.  Two  such 
plates  were  built,  and  this  procedure  was  used  for  configurations  X  and 
Z.  Temperatures  were  recorded  at  the  heating  face,  the  cooling  face, 
and  the  three  inside  stations. 

For  configuration  Y,  three  thermocouples  were  wired  in  parallel 
forming  a  multi- junction  thermocouple,  and  then  attached  to  an  aluminum 
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place.  The  aluminum  place  was  6  x  6  x  .03  inches  wich  chree  parallel 
sloes  cuC  chrough  ic.  One  sloe  was  3  inches  long,  and  bisecced  Che  6  x 
6  aluminum  place.  The  ocher  Cwo  sloCs  were  offsec  0.5  inches  on  eicher 
side  of  Chis  first  sloe.  The  chermocouples  lead  wires  were  placed  in 
these  sloes,  below  che  surface  of  Che  place.  Three  such  plates  were 
made  and  placed  at  0.5,  1.0,  and  1.5  Inches  from  che  heating  plate.  The 
temperatures  were  then  recorded  at  these  three  locations  and  at  the  heat¬ 
ing  and  cooling  face. 

Graphite-Epoxy  Composite 

The  graphite-epoxy  composite  samples  were  obtained  from:  Flight 
Dynamics  Laboratory  (Mr.  Robert  Achard) ,  Wright-Patterson  Air  Force  Base, 
Ohio. 

The  material  used  in  the  composites  was  manufactured  by  Hercules 
Incorporated,  Utah.  The  physical  properties  of  the  composite  as  supplied 
are  listed  below  (per  specification). 

Materials:  graphite  fiber,  epoxy  matrix;  3501-5A,  AS-1 

12  inch  prepreg  tape 
40  to  42  percent  resin  content 
fiber  density  (lb/in^)  0.0660 
resin  density  (lb/in^)  0.0460 
Composite  Properties: 

packing  arrangement,  random* 


*Data  were  determined  in  this  investigation.  All  other  values  on  this 
page  were  provided  by  the  manufacturer. 
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volume  ratio,  60+4  percent* 

62+3  percent 


fiber  orientation  0  degrees 

density  (lb/in3)  0.5785* 

(g/ca3)  1.604* 

void  content,  volume  Z  0 
Geometry:  length  (in)  6.0 

width  (in)  6.0 

thickness  (in)  0.5 

construction  105  lamina 


Temperature  Controller 

Manufacturer 

Model 

Type 

Range 

Supply 


Chrotnalox,  Lavergne,  Tennessee 
3804  PID 

J,  Iron- cons tan tan 

0  to  400°F 

120  VAC,  50/60  Hz. 


Data  Logger 

Manufacturer 

Model 

Type 

Scanners 

Scan  speed 

Channel  capacity 


John  Fluke  Mfg.  Co. 
2240A 

J,  K,  T,  or  R 
10  channels /second 
2.5  channels /second 
60 


*Data  were  determined  in  this  investigation, 
page  were  provided  by  the  manufacturer. 


Mount lake  Terrace ,  WA. 


All  other  vali as  on  this 
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Supply 


115/230  VAC,  60  Hz 


Digital  Volt  Meter 

Manufacturer 

Model 

Supply 

Matt  Meter  Transducer 

Manufacturer 

Model 

Type 

Supply 

Rating 

Thermocouple  Indicato_ 

Manufacturer 

Model 

Type 

Range 

Channel  Capacity 

Scan 

Supply 

Clock 

Manufacturer 

Model 


John  Fluke  Mfg.  Co. 
800 0A 

115  VAC/60  Hz. 


F.W.  Bell,  Columbus,  Ohio 
PT  2501F 
3  Phase,  4-Wire 
120  VAC,  50/60  Hz. 

.60  KVA 


Doric  Scientific 

D5-350 

J 

0  to  1600°F 
12 

Manual 

12  VAC,  60  Hz. 


Dime o-Gray  Co. ,  Dayton,  Ohio 
165 


Supply 


125  VAC,  60  Hz 


Variable  Transformer 


Manufacturer 

Type 

Supply 

Rating 

Thermocouple  Wire 

Manufacturer 

Type 

Diameter 

Description 

Fiberglass  Insulation 

Supplier 

Brand 

Dimension 

Thermal  Conductivity 

Heating  Elements 

Manufacturer 

Dimension 

Rating 

Brand 

Operating  Voltage 
Operating  Timetable 

Thermal  Grease 
Manufacturer 


The  Superior  Electric  Co.,  Bristol,  Connecticut 
216B 

120/240  VAC/60  Hz. 

0.98  KVA 


Leeds  &  Northrup,  Philadelphia,  Pennsylvania 
J;  Iron-cons  tan tan 
30  gauge 

Enamel  coating,  cloth  Insulation 


Service  Products,  Inc.,  Columbus,  Ohio 
#703 

2'  x  4'  x  2" 

.27  BTU-in/f t-hr-°F 


Electrofilm,  Inc.,  North  Hollywood,  California 

6”  x  6"  x  .007" 

2 

10  watts/in 
KAPINEX 

0  to  460  AC  or  DC 
-80°F  to  500° F 


Omega  Engineering,  Inc.,  Stamford,  Connecticut 
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Brand 


0T-201 


Maximum  Continuous 
Temperature 

Thermal  Conductivity 


500°F  (260°C) 

9.6  BTU-in/hr-f t 


Cooling  Plate 


Manufacturer 


Dimension 


In  Shop 
6"  x  6"  x  1" 


Material 


Aluminum 
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